rm(list=ls())

# Figure A.2

exp_claassen_m5 <- readRDS("claassen_mod5_output.rds")

pdf("fig_a2.pdf", height=4, width=4)

stan_rhat(exp_claassen_m5)

dev.off()
